Space-time versus particle-hole symmetry of collisions 
in a non-local kinetic equation for dense Fermi systems 
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Binary collisions in Fermi systems obey two fundamental symmetries corresponding to the space 
and time inversion and to the interchange of particles and holes. We show that beyond the local 
and instant approximation of scattering-in and -out integrals of a kinetic equation, only one of 
symmetries can be explicit while the other has to be covered by a constraint. This constraint, 
derived from the optical theorem, allows one to convert at need an explicit particle-hole symmetric 
form to the space-time symmetric form and vice versa. We implement this constraint to heavy ion 
reactions, where simulation algorithms require the space-time symmetry while former theories offer 
kinetic equations with the explicit particle-hole symmetry. 
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Realistic simulations of heavy ion reactions in the non- 
relativistic energy range have achieved a level on which 
basic problems with a reliable treatment of the Boltz- 
mann equation (BE) [|l) or of the quantum molecular dy- 
namics (QMD) |2J are already settled. The field is now 
open for a refinement of underlying equations of motion. 
Three kinds of improvements are at hand: momentum- 
dependent mean field B, medium effects on binary col- 
lisions ||], and non-local binary collisions (£J-0|. If im- 
plemented separately, each kind is covered by a corre- 
sponding theory. The mean field follows from the Landau 
theory of Fermi liquids [§||. The medium effects have 
been derived within the many-body statistics M]. The 
non-local collisions are supported by the theory of gases 
p[ . In realistic studies all the three kinds should appear 
together, but none of the mentioned theories covers this 
general case. 

The kinetic equation with the Landau-type mean-field 
and in- medium non-local collisions has been derived from 
non-equilibrium Green's functions [pL2|| . Although this 
kinetic equation includes all the three kinds of correc- 
tions, it is not suited for the numerical simulations be- 
cause of the symmetry reason. In the BE or QMD sim- 
ulations each single collision event plays a dual role of 
the scattering-in and -out so that the -in process natu- 
rally emerges as the space-time mirror of the -out pro- 
cess. Conversely, only theories offering non-local colli- 
sions with the space-time symmetry can be implemented. 
In the kinetic equation of |12[ the non-local corrections 
to scattering-out are not the space-time mirror of the 
scattering-in as a consequence of the particle-hole sym- 
metry of non-equilibrium Green's functions. Indeed, the 
particle-hole symmetry demands that hole-hole collisions 
(scattering-out) have the same scattering rate, includ- 
ing non-local corrections, as the particle-particle colli- 
sions (scattering-in). Accordingly, an explicit particle- 
hole symmetry excludes an explicit space-time symmetry. 

The explicit space-time symmetry is absent not only 



in fll2| but in any theoretical treatment of non-local scat- 
tering integrals [|13|-|l6|. The problem has not been rec- 
ognized, however, due to non-transparent forms of non- 
local corrections (in fact not aimed for simulations). We 
anticipate that in all mentioned approaches the collisions 
obey both symmetries, the particle-hole symmetry is ex- 
plicit while the space-time symmetry is hidden in the 
optical theorem. In this letter we derive a constraint 
with which one can interchange roles of symmetries mak- 
ing the space-time one explicit and the particle-hole one 
hidden. Although we focus on simulations of heavy ion 
reactions, the discussed problem exceeds merits of the 
nuclear physics as the symmetry of non-local collisions 
in Fermi systems is a general question. The key step of 
our approach can be easily translated to any formalism 
favoured in other fields of physics because it is a straight- 
forward substitution based on the optical theorem. 

To demonstrate the conflict between the particle-hole 
and the space-time symmetry, we discuss the system of 
hard spheres used in model studies of heavy ion reactions 
without H and with |3]|7) the Pauli blocking. Without 
the Pauli blocking the kinetic equation proposed by En- 
skog has non-local scattering integrals in which two col- 
liding particles are displaced by their diameter D, nw 



#4f^/-^-« + ] 



(i) 



Here / x = f(k,r,t), /+ = /(p,r + A 2 ,t),/ 3 = f(k - 
q, r, t),f 4 = /(p + g, r — A 2 , t), dP is the differential cross 
section, and A 2 = A--D is the displacement. The essen- 
tial part of the Enskog approach is a symmetry between 
the scattering-in and -out. If the scattering-out describes 
a process k,p —> k — q,p + q, its conjugate scattering-in 
process, k — q,p+q — ► k,p is its space and time mirror. 
Both symmetry operations are necessary. The time in- 
version reverses the process but it also flips directions 
of momenta which are then flipped back by space inver- 
sion. Due to the space inversion, the displacements in 



the scattering-in and -out integrals have opposite signs. 
Let us try to introduce the Pauli blocking of final states 
by hole distributions. From phase-space trajectories one 
expects that the final states have the same displacements 
as the initial ones, therefore the right hand side of (fil) 
should be extended as 



JdP [(l-/ 1 )(l-/ 2 -)/ 3 / 4 - - / 1 /+(1-/ 3 )(1-/+)] 



(2) 



This scattering integral violates the particle-hole sym- 
metry. Indeed, with the interchange of particles and 
holes the scattering-in and -out interchange their roles, 
therefore the directions of displacements appear reversed. 
This demonstrates the conflict between space-time and 
particle-hole symmetries. 

The conflict follows from incompatible concepts of the 
scattering-out in the classical and the quantum statis- 
tics. One has to recognize that a realistic collision has 
a finite duration and to compare these two concepts in 
the time picture. Within the classical picture of the 
scattering-out, the instant of the transition is attributed 
to the beginning of the collision when the particle en- 
ters particle-particle correlations. Related to the instant 
of transition, the collision process thus happens in fu- 
ture. A different picture is found within the quantum 
statistical approaches where the causal perturbative ex- 
pansion excludes any dependence on processes which are 
to happen in future. The scattering-out is hence rep- 
resented by a collision of two holes in which one of the 
holes annihilates the scattered-out particle. Since the 
particle is annihilated at the end of the hole-hole correla- 
tion, the collision process happens in the past, i.e., before 
the instant of the transition. The classical particle faces 
its partner of the scattering-out in front of it, while the 
quantum hole leaves its hole-partner behind. 

To break the ties of causality, we need the anti-causal 
expansion. The link between the causal and the anti- 
causal pictures of the binary collision is covered by the 
optical theorem for the two-particle T-matrix. We thus 
use the optical theorem to derive a constraint under 
which both pictures of the collision are equivalent. The 
implementation of the optical theorem can be done on a 
rather general level. The scattering-in and -out integrals 
result from anticommutators, {., .}, of the Kadanoff and 
Baym equation |||| 

{G > ,S<}-{G < ,S > } = {G > ,G > oT h (G < G < )T4 

~{G < ,G < oT R (G > G > )T A }. (3) 

Here, G < and G > are particle and hole correlation func- 
tions, the o denotes that G >:< closes one loop of the two- 
particle function on its right hand side. The T-matrices 
and pairs of single-particle correlation functions, (GG), 
obey standard two-particle operator products. With the 
first and second terms interpreted as the scattering-in 



and -out, respectively, the particle- hole symmetry is ap- 
parent from the interchange >< — ►<. 

Of the essential importance is that operators in (g) 
do not commute, what reflects the non-locality of colli- 
sions. The quasiclassical approximation of (0) thus yields 
non-local scattering integrals |12J. Naturally, non-local 
corrections to scattering integrals are determined by the 
order of operators. Formula (||) has the causal structure 
expressed by the order 'retarded-correlation-advanced' 
of the two-particle functions. Our aim is to rearrange (0) 
so that it will include the scattering-out in terms of the 
anti-causal expansion. In the optical theorem,!*] 



T R AT A = T A AT R , 



(4) 



on the left hand side, the time cuts of the retarded and 
the advanced functions restrict internal time integrals to 
the past while on the right hand side to the future. As 
we will see, the change in the order of operators results 
in the space-time inversion of non-local corrections. 

Further progress depends on the approximation of the 
T-matrix. Let us first approximate the T-matrix by 
Bruckner's reaction matrix for which the two-particle 
spectral function, A = (G > G > ), allows only for unoc- 
cupied states. The optical theorem then yields 

{G>, £<}-{G<, £>} = {G>,G> o T R (G < G < )T A } 

-{G < ,G < oT A (G > G > )T B }. (5) 

The first and second terms correspond to the scattering- 
in and -out integrals, respectively. Expression (0) has the 
desired explicit space-time symmetry contrary to explicit 
particle- hole symmetry of (0). 

To see it in detail, we substitute (^) or (||) into the 
Kadanoff and Baym equation and apply the quasiclassi- 
cal and quasiparticle approximations Mm- Keeping in- 
ternal gradients of the scattering integrals to the linear 
order, one arrives at the kinetic equation 



dfx dei dfx 
dt dk dr 



dei dfx 
dr dk 



dP-fcfz(i-h)(i-fc) 



dp ± (l-g)(l-f±)hf 



(6) 



* The theorem represents two alternative expressions of the 
anti-hermitian part of the T-matrix, M — i(Tr — Ta)- From 
the ladder equation in the differential form, T RA = V — Gr,a, 
with Gr.a given by the time cut of the spectral function A = 
i(0r -9a), follows ifT^ 1 -TJ 1 ) = -A. Multiplying M by 
T' 1 one finds T^M = i- iT~ x T A . Substituting T^ 1 = 
T A + iA one finds the causal side of the optical theorem 
M = TrATa- The anti-causal side, M = T A AT a , follows 



from MT- 



iTaT r x and the same substitution. 



where superscripts + and — correspond to (g) and (j|), 
respectively, and denote signs of non-local corrections: 
h = f(k,r,t), ft = f(p,r±A 2 ,t), ft = f(k±q±A K ,r± 
A 3 ,t±A t ), and ft = f(p+q±A K ,r±A 4 ,t±A t ). The 
dP + is obtained from dP~ by the flip of signs of all A's 
involved. More details about the limiting procedure and 
the differential cross section dP~ the reader can find in 
Jl2[ . One can see that with superscript — the scattering- 
out is the particle-hole mirror of the scattering-in, while 
with + it is the space-time mirror. 

All non-local corrections are given by derivatives 
of the scattering phase shift (p — lm lnT^(il, k,p,q,t,r) 

noun, 



A, = 



an' 



A K = 



\dcj) 
2dr' 



dp dq dk 



d dk' 
d(j> d<j) 
dk dq 



(7) 



The collision is of finite duration A t during which parti- 
cles can gain momentum Ax due to the medium effect 
on the collision. Three displacements A 2i 3,4 correspond 
to initial and final positions of two colliding particles, one 
position being fixed at r. The kinetic equation is justified 
to the linear terms in the A's which directly corresponds 
to the linear approximation in gradients. 

Without many algebraic details it can be indicated 
why the change of the causal picture into the anti-causal 
one results in the flipped signs of all non-local correc- 
tions. Writing the T-matrices as products of the ampli- 
tude and the phase, T R = |T|e l * and T A = |T|e _i *, one 
can see that the interchange of retarded and advanced T- 
matrices merely flips the sign of the phase shift <f>. As all 
A's depend linearly on 0, the non-local corrections to the 
anti-causal scattering-out have signs reversed compared 
to the causal one. 

Equivalency of both forms of kinetic equation (0) can 
be expressed by a simple constraint. Expanding the 
scattering-out to the linear terms, 

fdP ± (I- f±) (I- ft) Aft = lout ± J2 C outA Q , (8) 
■* a 

where I ou t is the local and instant scattering-out of the 
BE and C ou t are coefficients of non-local corrections, one 
finds that amplitudes of anti-causal and causal coeffi- 
cients equal while signs are opposite. Since both forms 
are equivalent, the sum of non-local corrections to the 
scattering-out vanishes, 






C« t A a = 0. 



(9) 



Perhaps, it is not necessary to remind that constraint 
@ is a direct consequence of the optical theorem (|J). 
The constraint gives us a freedom to express the non- 
local scattering-out as the particle-hole or the space-time 



mirror of the scattering-in, or if of interest, to suppress 
the non-local corrections to the scattering-out at all. 

Now the theory is ready for implementation into the 
BE or QMD simulations. In Fig. [TJ we present the re- 
sults of the BE simulation obtained with a modified BUU 
code of M . The non-local corrections are evaluated from 
formulas (R) and the Paris potential with no adjustable 
parameter. In accordance with (O), for the scattering- 
out we have used the standard local approximation of 
the cross section and the local Pauli blocking. In Fig. p] 
one can see that the production of high energetic pro- 
tons is enhanced due to non-local corrections. The effect 
of non-local collisions is more pronounced for larger im- 
pact parameter and for forward/backward angles where it 
achieves experimentally detectable values. A comparison 
with experiment and more technical details will be pub- 
lished elsewhere |17) together with the non-local QMD 
simulations. 
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FIG. 1. The angular resolved and absolute (lower right) 

BUU proton spectra for central collision of 209 Bi— > 136 Xe at 

55 MeV/A with and without non-local corrections. Three 

different impact parameter are chosen correspondingly. The 

local BUU code has been kindly provided us by W. Bauer. 

The conflict between the space-time and particle-hole 
symmetries deserves a discussion from a more general 
point of view. The above pragmatic solution aimed for 
simulations is not the full story because constraint (|9|) 
holds only for one particular case, the approximation of 
the T-matrix by Bruckner's reaction matrix. Moreover, 
when the particle-particle T-matrix T, evaluated from 
A = (G > G > ), is used in the scattering-in, the strict de- 
mand of the particle-hole symmetry implies to use hole- 
hole T-matrix T', evaluated from A' = (G < G < ) 1 for the 
scattering-out. This strict particle-hole symmetry cannot 
be demanded, however, within the Bruckner approxima- 
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tion because the conservation laws require equal scatter- 
ing rates of -in and -out processes but for \T\ ^ \T'\ the 
rates result different. 

The full particle-hole symmetry is obtained within the 
Galitskii-Feynman approximation pH] for which the two- 
particle spectral function A = (G > G > ) — (G < G < ) corre- 
sponds to the Pauli blocking of internal states (1 — /)(1 — 
/) — // = 1 — / — /. From the optical theorem (0) one 
then finds a modification of (@) 

fdp ± (i - ft ft)hft = h-f-f ± J2 c iV/ A - 

(10) 

with a constraint J^ a ^i-f-f^a = 0. This allows us to 
write two equivalent kinetic equations with either space- 
time or particle-hole symmetric scattering integrals 

d/i dei dfi _ ggi <9/i 
dt dk dr dr dk 

= /dP-/ 3 -/ 4 -(l-/i-/ 2 -) 

- fdp ± (i-ft-ft)hf£- (ii) 



Again, only one of the symmetries can be explicit, never- 
theless, the constraint guarantees that both symmetries 
are fulfilled at the same time. 

The kinetic equation (ttj\) is not suitable for numerical 
treatment by recent codes, because of negative scattering 
rates for 1 — / — / < 0. Recent codes deal exclusively 
with the Pauli blocking of form (1 — /)(1 — /) which 
keep the scattering rates positive. Numerical studies of 
non-local corrections are thus limited to the Bruckner 
approximation. 

We note that on the level of analytic theories, two 
equivalent kinetic equations (O) settle the conflict be- 
tween the space-time and particle-hole symmetries. The 
space-time symmetric form is more convenient for stud- 
ies of the conservation laws because it allows one to em- 
ploy symmetry operations developed within the classical 
theory of gases jy]. Since the non-local corrections to 
the scattering integrals make the conservation laws ap- 
preciably more complex, see |12|| , the explicit space-time 
symmetry offers a vital simplification. 

In summary, for non-local collisions the scattering in- 
tegrals cannot have explicit space-time and particle-hole 
symmetries at the same time. If one of the symmetries is 
made explicit, the other is covered by a constraint follow- 
ing from the optical theorem. Under the condition that 
the Pauli blocking of internal and final states are con- 
sistent, this constraint simply states that the non-local 
corrections to the scattering-out vanish. From the prac- 
tical point of view, this constraint gives us a freedom to 
select the explicit symmetry of non-local corrections in 
the scattering-out at convenience, or to suppress them. 



With recent simulation algorithms, tractable non-local 
corrections are restricted to Bruckner's approximation of 
the T-matrix. We do not find this limitation discourag- 
ing since Bruckner's theory has been quite successful in 
studies of the equilibrium nuclear matter. 
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